Microscopic model, spin wave theory, and competing orders in double perovskites 
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We present a microscopic theory of carrier-induced ferrimagnetism in metallic double perovskite 
compounds such as S^FeMoOe and S^FeReOe which have recently attracted intense interest for 
their possible applications to magnetotransport devices. The theory is based on an effective "Kondo- 
like" Hamiltonian treated here within the large- S expansion. We find that depending on the value 
of the carrier density the ground state is either a ferrimagnet or a layered antiferromagnet. The 
fcrrimagnetic state has a robust half-metallic electronic structure. The transition to antiferromag- 
netic phase is first order accompanied with the regime of phase separation. We study spin wave 
spectrum including quantum corrections and find strongly enhanced quantum effects in the vicinity 
of zero-temperature phase transition. 
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The recent discovery of the room temperature magne- 
toresistance in double perovskite compounds 1 has gener- 
ated intense interest in these materials because of their 
potential importance as magnetotransport devices as well 
as their rich and challenging properties. For applica- 
tions metallic fcrrimagnetic (FiM) compounds such as 
Sr2FeMoC>6 and Sr2FeRe06 are of particular interest, 
because both have high magnetic transition tempera- 
ture (T c ~ 420, and ~ 400 K, respectively) and their 
electronic structures have been suggested to be half- 
metallic. 1 Despite an enormous interest in half- metallic 
ferromagnets, most of the theoretical studies until now 
have been only by ab initio methods. 2-5 

In this letter we formulate a microscopic theory of the 
carrier-induced ferrimagnetism in this family of magnetic 
compounds. The theory is based on a minimal model 
Hamiltonian, in which only the lowest energy charge fluc- 
tuations, coupled to the local moments, are retained, and 
those which are moved to high energies due to strong cor- 
relations, are disallowed. The model has two parameters 
and allows for reliable analytical treatment of the physi- 
cally relevant situation. 

In the double perovskite structure Fe and Mo/Re ions 
form two face centered cubic sublattices and each pair of 
nearest-neighbor lattice sites is occupied by two different 
ions. In the ionic picture, Fe is in the 3+ valence state 
and half-filled d-shell forms a S = 5/2 local moment. The 
Mo is in 5+ state with one t^g electron in 4<i-shell. 2 This 
extra electron, which hybridizes with the same orbital 
states at the neighboring iron sites, is responsible for the 
system's lowest energy charge excitations. 6 In the classi- 
cal picture, only a spin down electron, with respect to a 
neighboring Fe local moment, can hop to this Fe site. For 
a spin up electron this hopping is blocked by the Pauli 
principle, as at the the iron site the spin up states are 
already all occupied. Thus, it is the exclusion principle, 
which couples the itinerant and the local spins kinemat- 
ically. This contrasts with the case of double-exchange 
systems, in which the intra-atomic Hund's rule is respon- 
sible for the coupling between the two subsystems. 7 In 
the magnetic double perovskites, as in double exchange 



magnets, ordering of local moments occurs due to the re- 
duction of the kinetic energy. The kinetic energy in the 
former case is minimized when the local moments are 
parallel to each other and antiparallcl to itinerant spins, 
and this causes the ferrimagnetism. 

Model Hamiltonian. - We now turn to a minimal 
model Hamiltonian reflecting the above discussion. We 
retain only low energy charge excitations in the sys- 
tem, which correspond to the processes (Fe d 5 ,B' d n ) <-> 
(Fe d 6 ,B' d 71 ' 1 ) (B' stands for Mo/Re). The Fe ions are 
treated as localized spins S = 5/2 and extra electrons 
from nonmagnetic Mo/Re ions as charge carriers. There 
are n = 1 and n — 2 carriers per unit cell in Mo and Re 
compounds, respectively. When a carrier is placed at a 
site with core spin S = 5/2, the total spin 5 can take two 
possible values 5 = 2 and 5 = 3. However, the maximal 
allowed spin for six electrons in a d-shell is 5 = 2. To 
project out the 5 = 3 spin state we introduce infinite 
local antiferromagnetic (AFM) coupling J — * oo between 
core and itinerant spins. 8 The carriers occupy three de- 
generate t2 g (d X y,dy Z ,dy Z ) orbital states. In the cubic 
lattice the ^-transfer matrix is diagonal in the orbital 
space and is non-zero only in the corresponding plane. A 
minimal Hamiltonian can be written as a sum of three 
two dimensional (2D) terms H = H xy + H xz + H yz , where 
each term corresponds to a given orbital state. All three 
terms have the same form given by (we further omit or- 
bital index for simplicity) 

H = -tj2 [4Jj° + H - c ] + J2 AcTn * - 5 Act ^' 

(ij)a " i j 

+ ^ J [S,Sj — An,] , (1) 

i 

where the first term describes an electron hopping be- 
tween nearest-neighbor Fe and Mo/Re ions, labeled by 
i and j, respectively. The operators di a (rii) and dj„ 
(hj) corrsepond to Fe and Mo/Re sublattices, respec- 
tively, 2A CT = E(Fed 6 , B'd°)-E(Fed 5 , B'd 1 ) is a charge 
transfer gap, and Sj (sj) stands for the localized (itiner- 
ant) spin. The last term with A = (S + l)/2 is chosen so 
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that the exchange energy vanishes for a state with 5 = 2. 

We now consider the ferromagnetic (FM) arrangement 
of local moments, which corresponds to the ground state 
of Hamiltonian (1) for carrier density n < n\ (see be- 
low). In the limit J — > oo , we restrict the Hilbcrt 
space to spin eigenstates with maximum allowed spin 5 
at a given site. Thus, the magnetic ion is in a state 
either with 5 = 5, without an extra electron, or of 
5 = 5 — 1/2, with an extra electron. The site states 
are labeled by |5, M) and M is z-projection of the spin 
5. Considering low lying excitations we can retain only 
the states with M = S and M = S — 1. Moerover, as 
the value of the core spin is relatively large S = 5/2 
we perform the 1/5 expansion around the ordered state. 
Then the local Hilbert space consists of four states: IS 1 , 5), 
\S,S - 1), \S - ±,5 - I), and \S - ±,5 - f). To 
describe transitions between these states we introduce 
magnon and fermion operators as B^\S, M) = \S,M — 1) 

and D]\S,M) = \S- 
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|), respectively. Through 
corresponding matrix elements, one can find the follow- 
ing relations d l = D l [l - B^B/AS - B^B/(32S 2 )] , and 
di = -D^yV^S to order l/S 2 . The state generated 
by D-[ corresponds to 5 = 5 + 1/2 and is projected out. 

As a next step, we express the Hamiltonian (1) in terms 
of the new operators, diagonalize its band part, and ar- 
rive in the momentum space to the following expression 
H = H + H 1 + H 2 , where 

Ha = ^2 { Ek [ a il ak l ~ & ki & k|J - Act^t^i} , 

k 

Hi = — == Yl {4- qT [^k, q a ki + M k , q fe ki ]B q + H.c.} , 

+ PkA^ + pgAM + p P% b li a Pi } . ( 2 ) 
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where = u k -D k j - v ^e lk!C d kl , fr ki = v k D kl + 

u k e ife *J ki , and u(v) k = ^[1 ± A C t/ E k }/2. The first 
term H , corresponding to the classical limit (5 — > oo), 
describes the band structure of the system. The elec- 
tronic structure of each orbital state is composed of the 
three bands: bonding (6 k j.) and antibonding (a k j.) of 
down spin electrons, and nonbonding (rf k f) of Mo/Re up 
spin electrons. _E k = \fs^ + A CT and £ k = 2t(cosk x + 
cos k y ) for the d^-orbital. For n carriers per unit cell 
each bonding band of a given orbital is n/3- filled and 
fully polarized. Given this band structure, the system 
is half-metallic for arbitrary model parameters (t, Act) 
and for any realistic value of carrier density. The last 
two terms Hi and H 2 [see Eq.(2)] are due to the quan- 
tum nature of the core spins. Spin up electrons, local- 
ized in the classical picture, can hop thanks to the quan- 
tum nature of the local moments, but leaving a trace 
of spin deviations along their trajectories. This pro- 
cess is described by H\. The last term H 2 describes 



fluctuations of the local moment generated by the hop- 
ping of spin down electrons. The corresponding ver- 
tices are given by A k!q = £ k - q u k , M kiq = £ k - q u k , 

P k°p = - £ kV k U p - £pU k U p , P^ b p = £ k M k W p + £pW k W p , 

and P£p = £ p u p u k - e k « k « P - 

Spin wave theory. - We now turn to analysis of the 
spin wave spectrum. The magnon Green function is ob- 
tained from the Dyson equation -D q , w = [u> — E q>w ] . 
Where magnon self-energy E = ^\ Ej and i = xy, xz, yz 
refers to a given orbital channel. We note that to the 
order we consider the magnon self-energy does not mixes 
the different orbital channels. This is due to the fact that 
the electron propagator as well as the magnon-electron 
vertecies are diagonal in the orbital space. To order 1/5, 
the self-energy E includes the second order perturbative 
contribution from Hi as well as first order contribution 
from H 2 (see first two diagrams in Fig 1.) (The analyti- 
cal form of E is lengthy and will be presented elsewhere.) 
Here we discuss the results. There are two branches of 
spin excitations, like those of localized ferrimagnets, the 
gapless Goldstone mode and the optical mode. 9 In addi- 
tion to this, there exists a Stoner continuum of spin-flip 
particle- hole excitations. 




FIG. 1. The leading order (1/5) and next-to- leading or- 



der (1/5 2 ) spin 



self-energies. Wave (solid) line stands 



for the magnon (electron) propagator. The small arrow de- 
notes an electron spin and gray (black) dot stands for the 
magnon-electron vertex proportional to \/\/~S (1/5). 

We below focus on the low energy mode. The spectrum 
of the magnons at the quasi-classical level is Heisenberg- 
like u q = 165Ji[l - 7i q ] + 165J 2 [1 - 72q], where 
Ji and J 2 are nearest and next-nearest neighbor ex- 
change couplings, respectively, 71 q = (cos q x cos q y + 
cos q x cos q z + cos q y cos q z )/3, and 7 2q = (cos 2 ^ + 
cos 2 q y + cos 2 q z )/3. The carrier-induced exchange en- 
ergies are given by J 1(2) = t 2 /(\QS 2 N) £ k 7i(2)™ k /^k, 
where 71 = 2 cos k x cos k y , 72 = cos 2k x + cos 2k y and n k 
is a Fermi distribution function. We also note that there 
is no Landau damping of spin waves at this order. 

In Fig. 2(a) we show carrier density dependence of 
the spin stiffness T> = 165(Ji + J 2 ) for different values 
of Act- At a fixed density the stiffness scales with the 
bandwidth: D ~ t for t » A CT and V ~ i 2 / A CT for 
t <C Act- As a function of density the spin stiffness ex- 
hibits strongly nonmonotonic behavior. It shows a max- 
imum at optimal filling n opt ~ 1 and then drops to zero 
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at some critical density n cr . 10 

We now analyze quantum corrections to the harmonic 
spectrum discussed above. These corrections are gener- 
ated by the magnon self-energy diagrams proportional to 
1/S 2 . 11 They are governed by fermionic excitations and 
are thus different from those in insulating magnets. At 
this order magnon self-energy includes fourth and second 
order contributions from Hi and H2, respectively, as well 
as mixed type terms. Graphically they are presented in 
Fig. 1. From the structure of the diagrams one easily 
verifies that renormalization of the electron propagator 
due to the scattering of spin waves as well as all relevant 
vertex corrections are consistently taken into account at 
this order. After collecting all the diagrams in this order 
one recovers the Goldstone theorem, that illustrates the 
consistency of the approximation. 




FIG. 3. Spin wave spectrum in [111] direction (in units of 
t = 1) versus momentum. Solid and dashed lines represent 
the harmonic and renormalized magnon dispersion law, re- 
spectively. 




Electron density 

FIG. 2. (a) Spin wave stiffness (in units of t = 1) versus 
carrier density, for different values of Act (Act =0, 4, and 
8, solid, dashed and dotted-dashed lines, respectively); (b) 
Ground state (GS) energy difference (in units of t, and for 
Act = 0) between AFM-II and ferrimagnetic states versus 
carrier density. 

The quantum corrections generate two main effects 
in the magnon spectrum: (i) they give rise to spin 
wave damping in the ground state, and (ii) they mod- 
ify the semi-classical dispersion law of magnons. The 
first nonzero contribution to magnon damping T q is pro- 
portional to 1/5 3 and is due to the spin wave scattering 
of carrier density fluctuations. To evaluate damping, one 
has to replace the bare magnon propagator in the self- 
energy diagrams by that obtained in harmonic approxi- 
mation (see Golosov in Ref. 11). In the long wave-length 
limit one finds T q cx (fk F }tq 6 /S 3 , where (/fe p ) denotes 
the Fermi surface average of /k = (sinfc 2 + sinfc^) 2 . 

Figure 3 shows magnon dispersion obtained semiclassi- 
cally (solid line) and with quantum corrections included 
(dashed line), for carrier densities n = 1, Fig. 3(a), and 
n = 1.8, Fig. 3(b). In the case of optimal carrier den- 
sity, n = 1, the quantum effects are small and practi- 
cally do not alter harmonic dispersion law. However, 
for a higher carrier density the quantum effects become 
strongly pronounced and result in a dramatic modifica- 
tion of the harmonic spectrum. They cause stiffening of 
the long-wave length excitations and strong softening of 
the zone boundary magnons [see Fig. 3(b)]. 



Ground state phase diagram. - As we have already 
mentioned, there is a critical carrier density n cr at which 
the spin wave stiffness vanishes [see Fig. 2(a)]. This indi- 
cates that there exist another type of magnetic ordering 
which competes with FiM order and becomes favorable as 
The most direct route to identify the symmetry 
of new magnetic order is to find the wave vector at which 
spin wave spectrum becomes unstable. The linear the- 
ory can not answer this question, because the harmonic 
dispersion vanishes identically in [111] direction, and no 
special soft mode is singled out at this order. The prob- 
lem can be resolved by considering quantum effects. The 
tendency of the zone boundary magnons to soften with 
increasing n [see Fig. 3(a) and 3(b)] can be interpreted 
as a precursor effect of a transition to a layered magnetic 
ordering with wave-vector Q = 7r/2[l,l,l]. The latter 
consists of alternating {111} FM planes. The magnetic 
moments of neighboring planes can be either misaligned 
by an angle 9, which is the case of canted magnetic struc- 
ture, or aligned antiparallel to form type-II antifcrromag- 
netic (AFM-II) ordering. Hence, there are two possible 
scenarios for the phase transition: (a) the transition is 
second order and is to a canted state or (b) transition is 
first order to AFM-II ordering. To decide which of these 
two is realized we have evaluated and compared the ki- 
netic energy of carriers in FM, canted and AFM-II back- 
ground of core spins. We have found that canted state 
never minimizes the energy and hence is never stabilized. 

In Fig. 4 (left panel) we show 2D (001) cut of AFM-II 
magnetic structure. In this magnetic structure carriers 
are confined within the ID FM stripes and have the dis- 
persion law ek — —2\/2t cos k (the factor \/2 is due to the 
two possible path of the charge transfer between nearest- 
neighbor magnetic ions). Figure 2(b) shows the difference 
between the energies (AE) of the AFM-II and FiM states 
versus carrier density for Act = (the same qualitative 
behavior is realized for finite values of Act)- As it is 
seen in this figure, with increasing carrier density AE 
shows the same trend, not coincidently, as the spin wave 
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stiffness [see Fig. 2(a)], and vanishes at carrier concentra- 
tion n cr slightly smaller than n cr suggested from linear 
spin wave theory. For n > h cr carriers gain more kinetic 
energy in AFM-II state. This consideration suggests the 
first order transition from fcrrimagnetic to AFM-II state. 
To investigate the possibility of phase separation accom- 
panying the first order transition we have evaluated the 
ground state thermodynamic potential ft = E — fin for 
homogenous ferrimagnetic and AFM-II phases as a func- 
tion of chemical potential \x. The homogenous state is 
destabilized with respect to phase separation when O 
of the fcrrimagnetic phase becomes larger than that of 
AFM-II state. The value of n at which this takes place 
gives the lower n\ and upper n 2 boundaries of carrier 
density for the phase separation regime. We find that 
the carrier density m, at which which the phase separa- 
tion occurs, is smaller then the critical carrier densities 
suggested form the above analysis, the following inequal- 
ity holds ni < h cr (n cr ) < n 2 . For A CT — we find 
m ~ 1.5 and n 2 ~ 2. 

o o. o o 
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FIG. 4. Left panel: Two dimensional (001) cut of AFM-II 
structure. Filled (opened) circles denote Fe (Mo/Re) ions and 
big arrows stand for the core spins. Black (grey) dashed lines 
indicate the allowed path for spin down (up) electrons; Right 
panel: ground state phase diagram of model (1). 

Based on the above analysis, in Fig. 4 (right panel) 
we present ground state phase diagram of Hamiltonian 
(1). For carrier density n < n\ the ground state is fer- 
rimagnct, for n\ < n < n 2 the systems phase separates 
(PS) in ferrimagnetic and AFM-II phases. The volume 
fraction occupied by FiM and AFM-II phases are given 
by vi = (n-ni)/(n 2 - ni) and v 2 = (n 2 - n)/(n 2 - m), 
respectively. Here, we point out, that in the real sys- 
tems, long-range Coulomb interaction, not considered 
here, most likely will prevent the phase separation at 
macroscopic scale resulting in the nanoscale domains. 12 
We have also analyzed the stability of homogenous AFM- 
II phase within the linear spin wave theory and found 
that this magnetic ordering is indeed stable for n> n 2 . 

So far, we have not considered the effect of direct 
AFM exchange J' between the local moments. At the 
semi-classical level, an effect of J' is to reduce carrier- 
induced FM exchange. However, in the compounds with 
large transfer gap Act , FM exchange can be largely sup- 
pressed and AFM exchange, if strong enough, can stabi- 
lize AFM-II magnetic ordering. Considering the metallic 
compounds with fully polarized carriers we have also ne- 
glected on-site Coulomb repulsion U between charge car- 
riers on 3d level of Fe-ions. (As at Mo/Re ions electrons 



reside in more extended Ad/bd level this interaction is 
much smaller and can be neglected. 2,4 ) This approxima- 
tion is fully justified for Sr 2 FcM06 compound that has 
one charge carrier per Fe-Mo cell and hence the Fe-ion 
is effectively quarter filled. Therefore the probability of 
generating doubly-occupied states is low, and the effect of 
the Gutzwiller projector on the uncorrelated wave func- 
tion is expected to be negligible. As for the Sr 2 FeReOg 
which has two charge carrier per formula unit, the on-site 
U, if it is not largely screened in the metallic state, may 
suppress the charge fluctuations and renormalize elec- 
tronic caracteristics of the system, such as an effective 
carrier density and their band-width. However, as far 
as magnetic properties of the systems are concerned, we 
expect that Coulomb repulsion can only change our quan- 
titative but not qualitative predictions. The spin wave 
spectrum as well as its carrier density dependence can 
be directly verified by inelastic neutron scattering exper- 
iments on Mo and Re based compounds. 
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